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Abstract. We show how to extend (and with what limitations) Avila's global 
theory of analytic SL(2,C) cocycles to families of cocycles with singularities. 
This allows us to develop a strategy to determine the Lyapunov exponent for 
extended Harper's model, for all values of parameters and all irrational frequen- 
cies. In particular, this includes the self-dual regime for which even heuristic 
results did not previously exist in physics literature. The extension of Avila's 
global theory is also shown to imply continuous behavior of the LE on the 
space of analytic M2(C)-cocycles. This includes rational approximation of the 
frequency, which so far has not been available. 



1. Introduction 

For an irrational /3, consider the quasi-periodic Jacobi operators on Z 2 (Z) 

(1.1) {H e ,^) k := v(9 + (3k)^k + c(9 + (3k)^k+i + c(9 + (3{k - l))^ . 

indexed by 9 G [0, 1). In this article, c and v are assumed to be functions on the 
torus, T := K/Z, with analytic extension through a band |Im(z)| < 5. As usual, 
v is real-valued which makes Hg-p a bounded self adjoint operator. We assume 



Operators of the form (1.1) arise in a tight-binding description of a crystal 
layer subject to an external magnetic field of flux /3 perpendicular to the lattice 
plane [T^l US]- In this context the functions c and v reflect the lattice geometry 
as well as interactions between the nuclei in the crystal; 9 represents a (random) 
quasi-moment um . 

A prominent example and the main motivation for this paper is extended 
Harper's model. Here, the electron is allowed to hop between nearest and next 
nearest neighboring lattice sites expressed through the couplings A2 and Ai,A3, 
respectively, 

(1.2) c x (x) := X 3 e- 27Ti(x+ ^ + A 2 + Aie 2 ™ (x+ * } , v(x) := 2 cos(2ttx) . 

We set A := (Ai, A 2 , A3) to simplify notation. The model is illustrated in Fig. [Tj 
Without loss of generality one may assume < A 2 , < Ai + A 3 and at least one 
of Ai, A 2 , A 3 to be positive. 
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Figure 1 . Extended Harper's model takes into account both near- 
est (A2) and next nearest neighbor interaction (Ai, A3). Relevant 
special cases are Harper's model (Ai = A3 = 0, also known as "al- 
most Mathieu") and the triangular lattice (one of A 1; A 3 equals 
zero). Extended Harper's model as given here was obtained by 
renormalization, reducing the number of coupling parameters from 
four to three. 



The model is sufficiently general to incorporate both rectangular and triangular 
(one of Ai, A 3 zero) lattice geometries. We also mention that Ai = A 3 = produces 
the well known almost Mathieu operator (in physics literature also known as 
Harper's model) with coupling X^ 1 . Contrary to Harper's model, not much is 



known about its extension given in (1.2) 



Analyzing the solutions to the generalized eigenvalue problem Hg.pip = Eip from 
a dynamical systems point of view leads to consideration of the analytic cocycle 
((3,A E (x)) where 



1.3) 



A E {x) :-- 



E — v(x) 
c(x) 



-ax 



0) 







We will call (/3, A E ) a J acobi- cocycle, generalizing the notion of Schrodinger cocyle 
associated with the special case c = 1. A precise definition of (analytic) cocycles 
is given in Definition |2.1 



A complication not present in the study of Schrodinger cocycles is the possibility 
of singularities. A Jacobi cocycle is called singular if for some xq G T, det A(xq) = 
c(x — (3)c(x) = (xq correspondingly is termed a singularity of (f3, A E )). 



By Oseledets' theorem Q22J; see also Appendix G), the asymptotic behavior 
of solutions to Hg.^ip = Etp is characterized by the Lyapunov exponent (LE) 
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Figure 2. Partitioning of the parameter space into regions I - III. 
So far duality based methods precluded a study of the self dual 
region III. 



associated with (/3,A E ). In particular, knowledge of the LE is a crucial tool to 
tackle the spectral analysis for the ergodic operators Hg.p. 

Our main result resolves the open problem of obtaining a complete description 
of the LE of extended Harper's equation as a function of the coupling A. Ear- 
lier attempts to compute the LE, both heuristic and rigorous, were based on an 
underlying symmetry of the Hamiltonian, known as duality [20[ [HI [HI |26l l2Tj . 
Relying on this symmetry, however, precludes the analysis of a certain significant 
region of couplings. 

With respect to duality, the parameter space splits into (see also Fig. [I]): 

region I: < A x + A 3 < 1, < A 2 < 1 , 
region II: < Ai + A 3 < A 2 , 1 < A 2 , 
region III: max{l, A 2 } < Ai + A3 . 

According to the action of the duality transformation, regions I and II are dual re- 
gions whereas region III is self-dual. For a precise meaning of duality for extended 
Harper's model we refer the reader to Appendix [F} A more general perspective 
on duality is given in J2U] • We reiterate that the duality based approach a priori 
excludes the self dual region. 

So far the only rigorous results on extended Harper's equation have been ob- 
tained in the dual regime (regions I & II) for /3 Diophantine (see (2.26) for a 
definition) [201 [El EE] . 

Ref. (21] quantifies the LE in region I for Diophantine /3. This formula for the 



LE in region I was obtained based on [16] which proves (spectral) localization for 
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Diophantine (3 in the interior of region I. Consequently, as shown in [16J, duality 
forces the LE in the interior of the dual region II to be zero. 

We mention that these results for the dual regime confirm non-rigorous com- 
putations by Thouless [TJIE6]. Since Thouless' considerations were also based on 
duality, they could not provide insight in region III either. 

As a consequence of their heuristic analysis of the Lyapunov exponent, in [TJ] 
Thouless and Hahn concluded localization within region I, and extended states for 
region II. Note that regions I and II can be viewed as extensions of the correspond- 
ing (metallic and insulator) regions of the almost Mathieu operator (Ai = A3 = 0). 

The self-dual region III, however, where next-nearest-neighbor interaction dom- 
inates, does not have analogues in the almost Mathieu family. In that region, the 
authors of [H] had to rely on numerical studies carried out for the case Ai = A3, 
which indicated singular continuous spectrum [14, 26J. We note that Ai = A3 cor- 
responds to the most physically relevant case of isotropic next-nearest-neighbor 
interaction, yet it turns out to be the most difficult one to study rigorously, because 
of the position of singularities of the associated Jacobi cocycle (see Observation 



3.1) 



The main achievement of the present article is to overcome the problem of 
analyzing the self-dual region of extended Harper's model. We prove the following: 

Theorem 1.1. Fix an irrational frequency (3. Then the Lyapunov exponent on 
the spectrum is zero within both region II and III. In region I it is given by the 
formula, 



;i-4) 



log 



log 



log 



'1 + Vl-4AiA; 

2Ai 



"1 + Vl-4AiA; 

2A 3 
l + v^ 



4AiA s 



A 2 + y/Xl - 4AiA 3 



ifM > A 3 , A 2 < A 3 + Ai 
if A 3 > Ai, A 2 < A 3 + Ai 
, if A 2 > A 3 + Ai . 



We emphasize that Theorem 1.1 holds for all irrational frequencies. 
Our strategy uses complexification of the analytic cocycle (/3, A E ) introduced 
in (1.3). For all A, where ((3,A E ) does not possess singularities, our method is 



based on Avila's global theory of analytic one-frequency Schrodinger operators 
[2], more precisely its extension to the GL(2,C) case, that we develop in Sec. 2.1 



A complication, however, arises in extended Harper's equation due to zeros of 
det A E (x) relating to zeros of c\(x). Singularity of the cocycle (/3, A E ) constitutes 
an important difference to analytic Schrodinger cocycles. [j] 

An important ingredient that allows us to both overcome singularities of the 
cocycle as well as extend our results to all irrational /3 is continuity of the Lyapunov 



1 We note that the holomorphic extension A(x + it) of the transfer matrix to all of C will 
enc ounter zeros for some some non-zero e for any A with (Ai,A3) 7^ (0,0), (see Observation 
3.1 1, however, as we will argue, these are much easier to deal with. 
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exponent in the energy, A, and (5. As will be shown, this is only a special case of 
a more general continuity property of the Lyapunov exponent, valid for the class 
of analytic cocycles. 

To this end denote by C^CIT, M 2 (C)) the class of 1-periodic functions on the 
real line, with analytic extension to some band, \\mz\ < 5, attaining values in 
the complex 2x2 matrices, M2(C). The space C W (T, M 2 (C)) is topologized by an 



inductive limit topology described at the beginning of Sec. 2.1 



Theorem 1.2. Given (3 irrational, the Lyapunov-exponent, L(f3 + .,.) : T x 
C W (T, M 2 (C)) — > R is jointly continuous. 

Remark 1.3. (i) It is important that (3 is irrational, since at rational /3 the 
theorem is known to be false even for the SL(2, C)-case (see e.g. Remark 
5 in [2] for a counter-example), 
(ii) If D = 0, L(P, D) = -oo for any /? G R. 



Theorem |1.2| is false when reducing the degree of regularity to even smooth 
cocycles. Furthermore, let N := N U {0}. For k G N U {oo}, let C k (T, M 2 (C)) 
denote the class of 1-periodic M 2 (C)-valued functions on the real line, topologized 
by the complete metric f\ 

/, eN /n m f max < i < fc max 3 . eT ||a- ? (D(x) - E(x))\\ , k G N , 

(1.5) pk(D, E) := < i Pn {D,E) ,_ nn 

yl^n=0 2" l+ Pn {D,E) ' h W ' 

where |.| is any fixed matrix norm . Here, k = is simply the continuous maps, 
formally d° = id, the identity operation on functions. For HN, d k D denotes 
the kth derivative taken componentwise. We restrict to the subset 

(1.6) C k := {D E C k (T,M 2 (C)) : log||D|| G L\T,dx)} , 



where the LE is finite (see (2.5)). By Fatou's Lemma C is open in C (T, M 2 (C)). 
Then, we can claim: 



Theorem 1.4. Theorem 1.2 is optimal, i.e. L{j3 + .,.): T x £°° — > R is discon- 
tinuous. Moreover, for any 6 No, + .,.): T x C w - s-IRzs discontinuous in 
the topology ofTxC k . 

Remark 1.5. (i) Certainly this implies discontinuity in C k (T, M 2 (C)), k G NU 
{oo}. 



(ii) Wang and You [2Z] have recently obtained Theorem |1.4 for SL(2, C)- 
cocycles by explicitly producing subtle counter-examples for any degree 
of regularity C k (T, SL(2, C)), k G N U {oo}. In case of M 2 (C)- cocycles, 
however, the construction can be done using elementary tools of harmonic 
analysis only, which is what is done here. 



2 For k €E No the metric is in fact derived from a norm, which however is not true for k = oo. 
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(iii) Discontinuity for k = is well known. In [12], Furman proved that for 
irrational /3, the LE is discontinuous at every (/3,D) where the limit (2.6) 
is not uniform. More recently, in [6], Bochi proved that cocycles with zero 
LE are dense in C°(T, SL(2,C)). 



Theorem 1.2 was preceded by a result in [18] which had already obtained con- 
tinuity of the LE for singular analytic cocycles but only under a Diophantine 
condition. The present formulation removes this limitation. 

The main accomplishment here, however, is to also prove continuity in the fre- 
quency at any irrational /3, which allows to consider rational approximates of the 
frequency. Such rational approximations are particularly useful for practical pur- 
poses since the resulting dynamical systems are periodic. For instance, continuity 
in this sense was the key assumption that Thouless used to compute the LE of 
extended Harper's model in region I. 

It is noteworthy that Theorem L2 arises as a consequence of ideas developed 
in Sec. [2] and illustrates the power of Avila's global theory approach [2J. The 
proof given here not only removes the problem of our earlier work to deal with 
Liouvillean frequencies but is also surprisingly simple. The key idea underlying 
this argument was communicated to us by Artur Avila (see more in the Acknowl- 
edgement). 

We organize the paper as follows. Sec. [2] forms the technical heart of our 
strategy to determine the Lyapunov exponent of extended Harper's equation. 
Based on analyticity of the cocycle (/3, A E ) we consider the complex extension 



A E (x + ie) =: A\ 



x , e G 



The underlying idea is to achieve " almost constant 



by considering the limits e — > ±oo. 

2.1, we extend Avila's global theory of analytic Schrodinger (or more 



cocycles 1 ' 
In Sec. 

generally, SL(2,C)) cocycles to the non-singular Jacobi (or GL(2,C)) case. In par- 
ticular, in Theorems 2.5, 2.9 and 2.10 we analyze the dependence of the Lyapunov 
exponent of the complexified cocycle on e, eventually enabling us to extrapolate 
to e = 0. 



As an immediate application, in Sec. 2J2 we establish continuity of the LE for 

This is the key ingredient to 



singular analytic cocycles as stated in Theorem 1.2 



deal with singularities of the cocycle. In the same section we also prove optimality 
of this result (Theorem 1.4). 

Having summarized some simple facts about extended Harper's model in Sec. 
[3j Section [4] obtains asymptotic expressions of its Lyapunov exponent. Sec. [5] 
carries out the extrapolation process leading to Theorem 1.1 In fact, as a result 



of our analysis we obtain the LE for all e (Theorem 5.1). 

An interesting property of the LE of extended Harper's equation is its indepen- 
dence of A2 in region III (see Fig. [T|. This symmetry is already revealed in the 
asymptotic analysis of Sec. [4] and hence inspires an alternative proof of Theorem 
1.1 1 which we give in Sec. [6j 
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We mention that the strategy we develop to obtain the LE for extended Harper's 
equation does not rely on the details of the functions c,v in (1.2); our method is 
hence applicable to a general singular quasi-periodic Jacobi operator with analytic 
coefficients. 

Finally, in view of future applications to other Jacobi operators, in Sec. [7] 
we point out differences and peculiarities of the Jacobi compared to 

Schrodinger cocycles which relate to possible zeros of det(A E (x)). As we will 
see, extended Harper's model provides useful examples and counterexamples. 

In particular, we address the question of almost reducibility which, for Schrodinger 
cocycles, is conjectured to provide a necessary criterion for (purely) absolutely con- 
tinuous spectrum in terms of the e dependence of L(f3,Af) near e = ("almost 
reducibility conjecture (ARC)" [21 13, EJ [1] ) . Extended Harper's equation provides 
a nice example that the ARC is false for the general Jacobi case due to possible 
zeros of det(A E (x)) on T. 

Acknowledgement. We are grateful to Artur Avila for his remarks on an 
earlier version of this paper where, among other things, he essentially provided 
a simple proof of continuity of the Lyapunov exponent of singular cocycles for 
all frequencies once such continuity for the Diophantine case is established, using 
the ideas of [2]. It turned out the same idea could be used to provide a simple 
proof of joint continuity (presented here), significantly simplifying our original 
approach. Additionally, our proof of continuity for the case of identically vanishing 
determinant follows his suggestions as well. We also thank Anton Gorodetski for 
useful discussions during the preparation of this manuscript. 



2. Complexified cocycles 
2.1. Jacobi cocycles - General results. We start by considering a general 



quasi periodic Jacobi operator with analytic coefficients as given in (1.1). We will 
be dealing with analytic functions on T. To this end, some topological remarks 
will come handy. 

Given a Banach space (X, \\.\\x), let Cf(T, X) be the space of X-valued functions 
on T with holomorphic extension to a neighborhood of |Imz| < 5, 5 > 0. We 
shall denote the band |Imz| < 5 about T by TV The set Cf(T,X) naturally 
becomes a Banach space in its own right when equipping it with the norm ||X||,j := 
su Pl ^ s \\X(z)\\ x ,ioTX(.)eC^T,X). 

To obtain statements independent of 5, we consider C W (T; X) := Lis>oCg(T, X). 
The appropriate topology for C W (T;X) is the inductive limit topology induced 
by \\.\\s. In this topology, convergence of a sequence X n — > X is equivalent to 
existence of some 5 > such that X n e Cf(T,X) holds eventually and \\X n — 
X\\s — > as n — > oo. 

If not specifically stated the sets Cg(T, X) and C W (T, X) will always be under- 
stood as topological spaces with respect to the above mentioned topologies. 
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For c(x) ytz o, we define the transfer matrix associated with the equation 
H e - ^ = Eip, 

{2 - 1} B {X} - c(x) { c(x) ) - c(x) A 

where solutions to the Schrodinger equation are obtained iteratively by 

B E > n {(3, 6) := B E [B + 0(n - 1)) . . . B E (6) , 

(2.2) B E ' n (j3, 6) := B E ' n ((3, 6 - n/3)" 1 , n > 1 . 

Equation (2.2) has a convenient dynamical formulation based on the following 
definition: 

Definition 2.1. (i) For (3 G K and D : T -> M 2 (C) Borel-measurable satis- 
fying 

(2.3) / log + ||L>(x)||dx < oo , 

we call the pair ((3,D(x)) a cocycle understood as a linear skew-product 
acting on T x C 2 defined by (x, t> ) i— > (x + (3, D(x)v). In (2.3), |.|| may be 
any norm on M 2 (C). If D G C a (T, M 2 (C)), a G N U {oo, w}, (/3, D(x)) is 
referred to as C Q -cocyle (also analytic cocycle if a — u). 
(ii) A C a -cocycle (j3, D(x)) is called singular if det -D(xo) = for some xo G T, 
in which case xq is referred to as singularity of the cocycle (f3,D(x)). 

Remark 2.2. For any D G C^(T, M 2 (C)) \ {0}, analyticity guarantees log||D|| G 
L 1 (T, dx) (for a simple argument see e.g. the proof of Lemma 2.9 in [18J). 

Since the functions c, v extend analytically to some band \lmz\ < 5, we may 
consider the complexified transfer matrix B E (x) := B E (x + ie) for |e| < 5. Set 

(2.4) 7 e (c) := / log|c(a; + ie)|dx = - / log |det Af (x)\dx , 



and write /o(c) =: 1(c) to simplify notation. 

Given a cocycle (/3, D) with (3 irrational, the Lyapunov exponent is defined by 



(2.5) 



L(/3,D) := lim - [ log|| J D (n) (x)||da; = inf - [ log||L> (n) ( 
n->oo n J j neN n Jj 



(2.6) = lim -log\\D^(x)\\ , 

n— >oo n 

where 

(2.7) D {n) (x) := D(x + (n - l)/3) . . . D(x) 

is the analogue of the n-step transfer matrix in ( 2.2[ ). Existence and a.e. indepen- 
dence of the limit in (2.6) follows by the sub-additive ergodic theorem. 
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In particular, for the cocycles ((3, A E ) and ((3,B E ), equation (2.1) implies the 
following relation 

(2.8) L(J3,B?) = L(I3,A?)-I e (\) , 

which allows us to focus on the analytic cocycle (f3, A E ). 

In view of Theorem L2, notic e th at by the sub-additive ergodic theorem the 
limit on the right hand side of (2.6) still exists a.e. (is L 1 and invariant under 
rotations by /3) even if /3 is rational but generally this limit will depend on x (and 
thus will not be equal to (2.5)). 



For rational ft, /3 — | with (p,q) = 1, we define the LE as given in (2.5) and 
note that by dominated convergence, 

(2.9) L(-,D) = - I logp (D (q \x)) dx . 

Here and later, p(A) denotes the spectral radius of a matrix A. 

Following, we denote the spectrum of Hg.p by £ and its a.s. components by 
S ac , and S sc , respectively. Our analysis of L((3,B E ), for E G E, is based on 
the complexification of the cocycle ((3, A E (x)). 

To this end, let ((3,D(x)) be a fixed analytic cocycle. We introduce the accel- 
eration, 

(2.10) */); f ):=iM« fl *')-" W '' 



2n h^o+ h 

The acceleration was first introduced in [2] for analytic SL(2, C)-cocycles. Ex- 
istence of u}(/3, D; e) is a consequence of convexity of L((3, D e ) w.r.t. e which is 
still true even if (/3, D) is singular. Albeit simple, convexity of L(f3, D € ) w.r.t. e 
for singular analytic cocycles will be shown to have far-reaching consequences in 
what is to come. We include a brief argument in Appendix [A] (see Proposition 



A.l). 



The next two theorems, Theorem 2J3 and 2^, extend statements proven by 
Avila for analytic cocyles with SL(2,C)-transfer matrices j2]. Avila's proof uses 
continuity of the Lyapunov exponent for Schrodinger cocycles established in [7] 
(see Theorem 2.15 below). 

It is straightforward to extend these statements to analytic cocycles whose 
determinants are bounded away from zero on a strip, making use of the following 
Lemma which we believe to be well known. For the reader's convenience, we 
provide proof in Appendix [E] 

Lemma 2.3. Let f e ^(TjC) with min xg Tr|/(x)| > 0. Then, there exists g G 
C"(R/2Z; C) such that g 2 = f . 

Remark 2.4. (i) From basic complex analysis it is clear that v7 can be defined 
holomorphically in a neighborhood of K, however it is not a priori obvious 
that this also yields a periodic function (which in general is 2- instead of 
1-periodic). 
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V 

(ii) The same proof shows that /«, can be defined in C w (IR/gZ; C). 



As a consequence of Lemma 2.3 Avila's results on analytic SL(2, C)-cocycles 
carry over to a non-singular analytic cocycle ((3, D) upon consideration of the 
"renormalized" SL(2, C)-cocycle, 

(2.11) D' :=—£=, D' eC u) (R/2Z;M 2 (C)) . 

In this context, it is useful to notice that a given matrix- valued function D : 
T M 2 (C) satisfying the hypotheses of Definition 2.1 (i) may also be considered 
as a function on K/2Z, in which case the Lyapunov exponents of the respective 
cocycles (for a fixed irrational /3) are related by a factor 2. 

Theorem 2.5 (Quantization of acceleration). Consider an analytic cocycle ({3, D(x)) 
where f3 is irrational and detD(x) bounded away from zero on a strip TV Then, 
the acceleration on Tg is integer-valued. 



Using Lemma 2.3 , we reduce to the original result for SL(2,C)-cocycles as stated 
in [2], noticing that 

(2.12) L{fi,D) = L(f3,D') + \ ! log|det(D)|dx . 



What is left to analyze is the second term on the right hand side of (2.12). We 
mention that this integral can be recast as the LE of a diagonal cocycle: 

(2.13) \ ^logldet^ld* = L (ft ( de 7' z > iet ° D(z) )) . 

In this case, however, the content of Theorem 2J3 may easily be checked directly, 
which is the subject of the following Lemma: 

Lemma 2.6. (i) For c £ Cf(T,C) with mm x£ T s \c(x)\ > the function I e (c) 
defined in (2.4) for e £ [—5, 6} is affine in e with derivative in 27rZ. 
(ii) If c £ Cf(T, C) \ {0} is not bounded away from zero, I 6 (c) is a piecewise 
linear, convex function in e with right derivatives in 2-k7L. 

Proof. (i) Since c(z) is holomorphic with no zeros on |Imz| < 5, log|c(z)| is 
harmonic on the same strip. Thus one computes, 

d 2 f f d 2 

— / log \c(x + ie)\dx = / — log \c(x + ie)\dx 

de J T J T oe 

f d 2 

(2.14) = - / — \og\c{x + te)\dx = , 



^ dx 2 

from which we conclude that I € (c) is affine on |e| < 5 as claimed. 

To show that the derivative is in 27rZ, we first consider the case when c is 
a trigonometric polynomial, i.e. c{z) = Ylk=-N Ckzk wnere z := e 2m ^ x+le \ 
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Since, 



(2.15) 



log|c(z) 



-N loglzl + lot 



2N 

E 

fc=0 



O-k-NZ 



(2.16) 



(2.17) 



it suffices to establish the claim for c(z) = Xlfclo c kZ k - 

Denote by Af(c; K) the number of zeros of c (counting multiplicity) on 
a compact subset K of C and let 3(c; K) denote the associated zero set. 
Since c is bounded from zero on T,5, we have 

AT (c; D(0,e~ 2 ^ = AT (c; D(0, e - 2 ^ , 

whenever < |e| < 5. 

Using Jensen's formula we thus conclude for < |e| < 5, 



\og\c(x + ie)\dx = -2neM(c; {0}) + V] n(z)\o£ 

ze3(c;D(0,e- 2 ' r ' s )\{0}) 



-2nt 



+ d 



where d = d(c) e C and n(z) is the multiplicity of z if z is a root of c and 
defined zero otherwise. 

For a general c, uniformly approximate c on T<j by trigonometric polyno- 
mials c n . As shown in [18] (see also Remark 2.16), I e (c n ) — >■ J e (c) uniformly 
on [—5, 5]. 

In summary we obtain lines, I e (c n ), with slopes in 27rZ converging uni- 
formly on [—5, 5} to the line I e {c), whence forcing the derivative of the limit 
to be in 27rZ (see also Fact |2.17[ ) . 
(ii) To prove statement (ii), factorize c on according to its roots 



(2.18) 



c(x) = h(x) Y[ (e 2nx - e 2 -(^+^)) nj , 
i=i 

where nj is the multiplicity of the jth root and h is zero-free and holomor- 
phic on TV 

Thus, making use of (i), 



(2.19) 
(2.20) 



He) 



n „ 

Vrij / \og\e 2ni{x+ie) - e 2ni ^ +ie ^ \ + 1(h) 

n „ 

Vn 3 / \og\e 2ni{x+ie) -e 2ni{x i +ie i ) \+2irNe 

.7=1 1 



for some iVGZ. 
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Applying Jensen's formula separately to each summand on the right 
hand side of (2.20), we conclude 

(2.21) / l g|e 27ri(x+ ^ - e 2ni ^ +ie ^\ = -27rmin{e,e j } . 

Combining (2.20) and (2.21) yields the claim. 

□ 

The following special case will turn out to be of relevance for our further devel- 
opment: 

Definition 2.7. [2J If L((3, D e ) is affine about e = 0, the associated cocycle (/3, D) 
is referred to as regular. 

Next we explore the relation between regularity and the dynamics induced by 
the cocycle (f3, D). To this end we define: 

Definition 2.8. An (analytic) SL(2,C)-cocycle (/3,D) is called uniformly hyper- 
bolic if there exist (analytic) maps s, u : T — > PC 2 such that 

(i) D(x)u(x) = u(x + ft), D(x)s(x) = s(x + 0), 

(ii) for w e C 2 , ||w|| = 1: ||-D(a;)iu|| > 1, if ir(w) = u(x), and ||Z?(x)iu|| < 1, if 
ir(w) = s(x), Wx G T. Here, n denotes the canonical projection of C 2 onto 
PC 2 . 



Using Lemma 2.3 , the following is obtained as a mere corollary of Theorem 6 
in 0: 

Theorem 2.9. Consider an analytic cocycle (/3, D) where /3 is irrational and 
det-D bounded away from zero on \lmz\ < 5. Assume L((3,D') > 0. Then the 
cocycle ((3, D') is regular if and only if it is uniformly hyperbolic. 



In particular, applying Theorem |2.9| to Jacobi operators, we obtain a statement 
which will be central to compute L(/3, B E ) on the spectrum: 

Theorem 2.10. Let (5 irrational and assume detA E (z) bounded away from zero 
for \lmz\ < 5. If E e £ with L(B E ,f3) > 0, the cocycles (/3, A E ) and (/3, B E ) 
cannot be regular. 

We note that 

(2.22) HJ3,{A?Y) = L(J3,B?) . 



The proof of Theorem 2.10 is a consequence of the following two Lemmas and 
Theorem 12.91 

Lemma 2.11. Let (3 irrational and assume det A E (z) bounded away from zero 
for |Imz| < 5. Then, (/3, (A E )') is regular if and only if (f3,A E ) is regular. 
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Proof. Lemma 2.6 implies that for |e| < 5 we have 

(2.23) L{fi, (Af )') = L(/3, Af) + 2nNe + T , 

some r e M and iVeZ. Since an analytic cocycle is not regular if and only if its 
acceleration has a jump discontinuity at e = 0, we obtain the claim. □ 

Lemma 2.12. Let (3 irrational and assume det A E (z) bounded away from zero 
for \lmz\ < 5. If E 6S, the cocycles ((3, (Af)') and ((3,B E ) cannot be uniformly 
hyperbolic. 

The statement is well known for (/3, B E ). The claim for (/3, (Af)') is obtained by 
similar means, in addition making use of ergodicity. We give a proof in Appendix 
O 



Theorems 2.5 and 2.10 form the core of our method to determine the Lyapunov 
exponent of extended Harper's model. They characterize the e-dependence of the 
Lyapunov exponent of the complexified cocycles ((3,Af) and ((3,Bf) for energies 
in the spectrum. 

2.2. Continuity of the LE for singular analytic cocycles. Before applying 
the results of the previous section to extended Harper's equation, we present a 
surprisingly simple proof of the continuous dependence of the LE on the cocycle 
upon variaton over the analytic category, further illustrating the use of complex- 
ified cocycles. 

First recall the following basic, however very useful fact, even valid for contin- 
uous M 2 (C)-cocycles: 

Theorem 2.13. L :T x C°(T, M 2 (C)) — > R U {— oo} is upper- semicontinuous. 

Remark 2.14. As an immediate corollary, we have that L : T x C°(T,M 2 (C)) — > 
]RU{— oo} is continuous at every D) with L((3, D) = — oo. By the same reason, 
L : T x C°(T,SL 2 (C)) ->• R 0+ is continuous at every (J3,D) with L(f3,D) = 0. 
Here, R 0+ :={i6»:i> 0}. 

Because of its usefulness and since the argument is very short, for the reader's 
convenience, we include a proof of Theorem 2.13 in Appendix [Dj 



Given 5 > 0, the following subsets of Cg(T, M 2 (C)) will be of interest in the 
subsequent discussion: 

(2.24) A.g(T, M 2 (C)) := {D G C%(T, M 2 (C)) : det D(x) £ 0} , 

(2.25) 13$ (T, M 2 (C)) := {D G C?(T, M 2 (C)) : A/"(det D; T) = 0} . 

We note that A£(T,M 2 (C)) is open in the Banach space C^(T,M 2 (C)). As 
before, it is useful to consider A"(T,M 2 (C)) := U S>0 A^(T, M 2 (C)) which is open 
in ^(T, M 2 (C)) relative to the inductive limit topology introduced at the be- 



ginning of Sec. 2.1 Analogously, one defines S W (T, M 2 (C)), which is open in 
A W (T,M 2 (C)). 
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In [TS], we proved that for a given Diophantine /3, L(f3,.) is continuous on 
Ag(T, M 2 (C)). We recall that is called Diophantine if there exists < b{0) and 
1 < r(/3) < +oo s.t. for all j G Z \ {0} 

(2.26) |sin(2vrj/3)| > ^ . 

The Diophantine condition was imposed in order to deal with singularities of the 
cocycle. Even though it was speculated in [18] (and in an earlier, related result 
[T7] ) that the theorem should hold true for all irrational (3, actual proof did not 
follow from the method of [18]. 

In particular, the theorem proven in [18] does not imply continuity in the fre- 
quency /?. For practical purposes, however, continuity upon rational approxima- 
tion is most desirable since the LE of rational approximates takes a simple form 



due to periodicity of the cocycle (see (2.9)). The main result of this section 



Theorem 1 1. 2 [ also implies continuity in this sense. 

Questions of continuity of the LE have been actively studied in the recent past; 
for a brief survey we refer the reader to the introduction of [T8] 127] . For the 
present development we will make use of the continuity statement obtained for 
non-singular cocycles as proven in |7j; for further use we state this result here: ^ 

Theorem 2.15 Q7J). The Lyapunov exponent L{(3 + ., .) : T x £ W (T, M 2 (C)) R 
is jointly continuous at every irrational /3. 



Let us postpone for a moment the proof of Theorem 1.2 and first show that 



continuity in the analytic category is the best one can expect (Theorem 1.4). 



Recalling (2.13), it is enough to proof the following statement: 

Proposition 2.1. (1) The function I : C°°(T, C) R U {-oo} ; defined by 

(2.27) 1(c) := [ log|c(x)|dx , c G C°°(T,C) , 
is discontinuous at 

(2.28) c(x) 



T 



e -l/y/\sm(^x)\ ;i f x> Q ) 

, otherwise. 



Here, C°°(T, C) is topologized in analogy to (1.5). 
(2) For k G N , / : C W (T, C) ->• R is discontinuous at c(x) = sin 2fc (27rx) in the 
topology ofC k (T,C)- 



Remark 2.16. (i) Based on Jensen's formula, we showed in [18] (see Lemma 
2.9 therein) that I is continuous when defined on C W (T, C) (in inductive 
limit topology). 



3 In [7], Theorem 2.15| was stated and proven for the Schrodinger, SL(2,M) case; strictly 
speaking, the extension to SL(2, C) follows from [17]. The obvious generalization to non-singular 
cocycles has been carried out explicitly in Sec. 4 of [27] . 
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(ii) Concavity of the log and bounded convergence implies that for any k G 
No U {oo}, I is continuous at any c bounded away from zero. 

Proof. Let (p(x) be the standard C°° bump-function with supp</> G [—2, 2] and 
(j)(x) = 1 for x G [-1, 1]. Set 4>{x) := 1 - <j>{x) and let a k = ||0 (fc) ||[_ 2 ,2]- 
Consider the sequence of C°° functions on T given by 

(2.29) - f n := c(x) [<p(4n(x - §)) + 0(4nx) + <p(4n(x - 1))] 



(2.30) 



+c( 



An ' 



4>(4n(x - §)) + 0(4nx) + <p(4n(x - 1)) 



Claim 2.1. Poo(/n,0) — )■ as n — >■ oo. 



Proof. Since the remaining terms in (2.29) can be treated similarly, we will focus 
on showing that 

(2.31) p oo ((j)(Anx)c(x),0) — > , as n — > oo . 



By induction, for fceNwe have for x^O 

i 

e V\* in2 ™\ | sin 27rx|~ 3fc p6fc-i( \f\ sin27ra;|), 



(2.32) 



c^(x) 



(2.33) c (2fc+1) (a;) = e | sin27ra;|" 3{/c+1/2) p 6fc (v / l sin27ra;|) cos2tt.t, 
and c^(0) = 0, where pk(x) is a polynomial of degree k. Thus Vd, fceN, 

(2.34) nic^z)!!^^] =: n d 6™ ->■ , as n ->• oo . 

Given e > 0, choose K G N such that X]fc>i^ ^ fc < e - The claim follows if we 
can argue that 

K 

(2.35) 



^2 2 -k_Pk((J)(4:nx)c(x),0) 

k=l ' 



K 



p k ((j)(Anx)c(x),Q) 



< 



J2^ k pk{<P{^nx)c{x),0) ^0 



n — >• oo . 



k=l 



Indeed, by (2.34) for each k G N we have 

k 

(2.36) ||9 fe (c(x)0(4nx)) || < ^ 

p=0 



a p (4n) p bl -> , n — >■ oo . 



□ 



Approximate /„ by trigonometric polynomials p n satisfying ||/ n — p n || fc < 



1 \n 



and set c n = c + p n . For n G N, let J n := [0, 4-) U [\ 
We estimate: 



j_ i 

4n' 2 



+ [1-^,1)- 



(2.37) 



I{c n ) < -log(2) + / log|c n |dx 



T\J„ 



for all n G N. Since the integral on the right hand side of (2.37) converges as 
n — > oo to 1(c) > — oo, we have I(c^) 1(c). 



The second statement follows in the same way with c(x) = sin (2nx). 



□ 
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Proof of Theorem 1.2 The key observation, based on the idea of Artur Avila, is 
that by convexity of the LE in e, it suffices to show continuity of the LE away 
from T, i.e. on < |e| < 5 for some small 5 > 0. 

To see that this indeed is sufficient, suppose for a moment the LE was continuous 
away from T. Then for any 5* with < 5* < 5, L(f3 + r, B € ) — > L(/3, A e ) uniformly 
on 5* < \e\ < 5 as r ->■ and \\B - A\\ s -> 0. 

Lemma 2.17. Consider J 7 := {/ : [0,1] — > R convex} as a closed subset of 
(C([0,1],R), ||.||[ ,i])- Forx E (0,1) denote by D±(/)(x) the right(+) and left(-) 
derivative of f E T, respectively. Given K C (0,1) compact, define T^(f) := 
sup xeK \D±(f)(x)\ for f E J 7 . Then, T% is locally bounded. 

Remark 2.18. This is just a more elaborate version of the following basic fact 
for a sequence /„ : (0, 1) — > R of convex functions: If f n — > f pointwise, then 
£>_(/) < liminf^ D4f n ) < hmsup^ D + (f n ) < D + (f) (e.g. [25]). We 



prove Lemma 2.17 in Appendix [E 



Thus, the acceleration is locally bounded, i.e. there exist 77, k, K > such that 
whenever \\B — A\\g < rj and |r| <Kwe have 

(2.38) sup \u(P + r,B;e)\ < K , 

5*<|e|<<5 

which since u monotonically increases in e yields 

(2.39) sup\u(/3 + r,B;e)\ < K . 

Continuity of the LE then follows upon successive approximation, i.e. 
\L(f3 + r,B)-L(f3,A)\ < \L(/3 + r, B) - L(f3 + r, B s *)\ 

+ \L{p + r, B s .) - L((3, A s *)\ + \L{p, A s *) - L{fl, A)\ 
(2.40) < 2KS* + \L(/3 + r, B 5 *) - L(/3, A s *)\ , 

for every \\B — A\\ < rj and |r| < K. 



Since 5* was arbitrarily, Theorem 1.2 would follow provided that the second 



term in (2.40) can be made arbitrarily small. 



For ^4^(T, M2(C)) this is an immediate consequence of Theorem 2.15, since for 
5 > sufficiently small, det A(z) + if < |Im(z)| < 5. Thus, 



Theorem 2.19. Given (3 irrational, Theorem 1.2 holds for the restriction of L((3+ 
.,.) to T x A 0J (T,M 2 (C)). 

What is left is to consider is the case when the determinant vanishes identically 
on T. From a dynamical point of view this situation is analogous to the case of 
uniformly hyperbolic cocycles (with obvious abuse of terminology since the deter- 
minant is identically zero), where it is known that the LE behaves continuosly. 

That this reasoning may be used to extend above results on continuity to the 
case of identically vanishing determinant, was also pointed out to us by Artur 
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Avila which we gratefully acknowledge. Following we give the details to this line 



of argument, which combined with Theorem 2.19 completes the proof of Theorem 

o 

Fix a cocycle {(3, A) with irrational and detA = on T. Clearly, this auto- 
matically implies det A = on T$, the domain of holomorphicity of A. 

Consider first the case of non-trivial A (i.e. A ^ on T). Continuity for trivial 
A is dealt with in Lemma [2.201 

Without loss of generality we may then assume that the first row of A does not 
vanish identically on T. Thus, there exists < 5* < 5 such that for all z £ Tg* \T, 
(An(z), Ai 2 {z)) 7^ 0, whence dim Ker.A(z) = dim Rany4.(z) = 1 on IV \ T. 

Letting, 

(2.41) u(z) := (-^»W) , .(,) := A(,) 

we thus obtain analytic functions defined on Tg* \ T with Spa.n{u(z)} = Ker A(z) 
and Span{t>(z)} = Ram4(z). 

For z e Tg* \ T, define C(z) := (v(z — j3)\u(z)). We distinguish the following 
two cases: 

Case 1: det C(z) ^ on TV \ T. Then, making 5* > sufficiently small, 
KerA(z) and Ram4(;? — (3) are transversal i.e. C(z) £ GL(2, C) for all 
< |Im^| < 5*. 

For z G TV \ T we thus obtain an analytic conjugacy of the cocycle 
(f3, A) given by 

(2.42) C(z + P)- x A{z)C(z) = ^ °) =: D(z) . 

Here, c(z) is holomorphic and non-zero on TV \ T. 
We thus obtain, 

(2.43) L(/3, A e ) = L(/3, D e ) = [ \og\c{x + ie)\dx , 



T 



for < |e| < 5*. From above construction, it is clear that c(z) depends 
continuously on both A and (3. 
Case 2: If det C(z) = on \ T, Ker^z) and Ram4(z — /?) coincide on 
IV \ T. In this case A(z + (3)A(z) = on T 5 *. Thus L((3, A e ) = -00, 



\e\ < 5* which directly leads to continuity by Theorem 2.13 
Since in Case 1 (c(z) 7^ on < \Imz\ < 5*), 

(2.44) L((3,A t )=L((3, ^ x + ^ % = jf i og | c ( x + ie)\dx + L(/3, (l ° 
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we conclude that for < |e| < 5*, stability of the LE under any analytic perturba- 

•»■-»-—■ »—(::: 

More generally, the following is true: 

Lemma 2.20. The Lyapunov exponent L(., .) : T x C°(T,M 2 (C)) — > 1R U {-oo} 
is (jointly) continuous at any constant cocycle. 

Remark 2.21. Note that similarly to the continuity of the LE over the uniformly 
hyperbolic cocycles, the frequency does not have to be irrational. 

Proof. Let (fl,A) be a constant cocycle. We distinguish three cases depending on 
the relative magnitude of the eigenvalues of Ai, A2 (counting multiplicities) of A. 

The case |Ai| 7^ |A 2 |, was examined (in greater generality) by Ruelle in [23] . 
We mention that Ruelle's formulation only states continuity (in fact even real- 
analyticity) in the matrix-valued function, however it is clear from his proof that 
the same strategy shows joint continuity in (j3,B). 

If Ai = A 2 = 0, L(f3,A) = —00 whence continuity follows trivially from upper- 



semicontinuity, Theorem |2. 13 

Finally, if |Ai| = |A 2 | 7^ 0, det A 7^ 0, the remark to Theorem |2.13| immediately 
implies continuity of the LE for the "renormalized" cocycles 
(/3, A/ y |det(A)|). This, in turn implies the claim by Remark 



{p,B/y/\detB\) at 
□ 



2.16 
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□ 



3. Extended Harper's model 

To simplify notation, throughout the paper we shall suppress some of the depen- 
dencies of Hxfi.fi if the context permits. Following we summarize some statements 
necessary to apply above results for general Jacobi cocycles to extended Harper's 
model. 

First we mention: 

Observation 3.1. Letting z = x + is, c(z) has at most two zeros on the strip 
< x < 1; the zeros occur if at least one of the following conditions is satisfied: 

(3.1) A ie - 2m = A 3 e 2m , 

(3.2) A ie - 27re + X 3 e 2ne = ±A 2 . 

Necessary conditions for real roots are X± = A3 or Ai + A3 = A 2 . Moreover, for 
Ai = A3, c(z) has real roots if and only if 2A3 > A 2 . If Xi ^ A3 real roots occur if 
and only if X\ + A 3 = A 2 . 
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Application of Jensen's formula yields (see [TB] for e = 0) 
(3.3) 

'log A 3 + 27T6 



log Ai — 27re 



log 
Jog A 2 



2A!A 3 



A2 + \/Ai~ 



4AiA 2 



if A 3 e 2 ^ > Aie" 27re > 
and Aie" 2 ^ + A 3 e 27re > A 2 > , 
if Aie~ 27re > X 3 e 27Te > 
and Aie" 2m + A 3 e 2?re > A 2 > , 

if Aie" 2m + A 3 e 2 ^ e < A 2 and A 1; A 3 ^ 

if Aie -2 ™ + A 3 e 2 ^ e < A 2 , Ai or A 3 = . 



where we set I(c\) =: J(A) in (2.4). 



We note that (3.3) explicitly shows that I e (X) is jointly continuous in e and A. 



Application of Theorem 1.2 to the model under consideration results in 



Theorem 3.1. For fixed irrational (5, the Lyapunov exponents L((3 + r; Bf , A) 
and L((3 + r; Af, A) are jointly continuous in e, X, E and r. 



We emphasize that Theorem 2.5 and 2.9 require analytic cocycles with deter- 



minant bounded away from zero. As discussed in Observation 3.1, for any A 
with Ai 7^ A 3 or A 2 7^ Ai + A 3 , c(z) is non-zero on the strips {|Im2:| < ex} and 
{|Imz| > e 2 }, for some < t\ < e 2 . 



Thus the hypotheses of Theorem 2.5 and 2.9 apply to (/3, Af) for small and 



large e (in magnitude) when excluding Ai = A 3 and A 2 = Ai + A 3 . Let 
(3.4) K := {A : X 1 ^ A 3 , A 2 ^ X, + A 3 } . 

Notice that upon exclusion of the two planes described by 7Z C , we do not loose any 
information; the derived statements about L(f3, B E ) will extend to all values of A 



using Theorem 3.1 as well as continuity of the spectrum in the Hausdorff metric. 



4. Asymptotic analysis 

In this section we aim to obtain an expression for L(f3,Af) valid for large |e| 
and any choice of A. We mention that this strategy can be employed generally to 
Jacobi cocyles consisting of trigonometric polynomials. 

The basic idea is to reduce the non-trivial problem of computing the LE of a 
given non-constant cocyle to an "almost constant" cocycle by taking |e| — > 00. 



By Theorem 1.2 the LE is stable under small analytic perturbations, whence the 
task boils down to consideration of a constant cocyle whose LE can be computed 



trivially. Finally, quantization of acceleration, convexity and the Theorems |2.5 
and 2.10 allow us extrapolate to e = 0. 

Our strategy to determine the Lyapunov exponent was motivated by consider- 
ations of Avila's in [2j. There, complexification is employed to prove the Aubry- 
Andre formula (see Appendix A in [2]). 
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Since the present model specializes to the almost Mathieu equation when Ai = 
A3 = 0, as a warmup, we present a proof of the Aubrey- Andre formula alternative 
to Avila's in that it replaces a geometric argument with one along the line of what 
is to come. 

Setting Ai = A3 = 0, we obtain the transfer matrix of the almost Mathieu 
equation 

E — 2/i cos(27rx) — 1 N 
1 



(4.1) 



B E (x) 



To simplify notation we put fj, — A 2 ■ Here A E = B E 
Cc 

(4.2) 



Complexifying (4.1), as e -> +00 
B E (x) 



^e~ 2wix e 2m 







: fie 



- 2nix e 27re DJx) 



-l + o(l) 
o(l) 

uniformly in x G T. Here, the the o(l)-terms are perturbations by trigonometric 
polynomials with only non-negative harmonics. We mention that for SL(2, R) 
cocycles the reflection principle implies that L(f3, B E ) is an even function in e. 
Applying Theorem 1.2 to (4.2) we conclude, 

(4.3) L(P,B?)=]og\fi,\ + 2n\e\, \e\ > e' , 

which gives the desired asymptotic formula for the complexified Lyapunov expo- 
nent of the almost Mathieu equation. 



Remark 4.1. Making use of Theorem 2J3 and 2.10 the asymptotic relation 
(4.3) can easily be extrapolated to e = producing the complex Aubry- 
Andre formula 

(4.4) L(p, B E ) > max{0, log \fx\ + 2vr|e|} , 

with equality holding for E e S. For details we refer to [2], Appendix A. 

Following we shall focus on the case when at least one of Ai, A3 is nonzero. 

Proposition 4.1. For at least one of X±, A3 nonzero, 

L(P,Af) = log 



1 + X /1-4A 1 A 3 



+ 2ir\e\ , e > e , 



(4.5) 
some e > 0. 

Proof. Case 1 - both A x and A3 are nonzero: Then, uniformly for x G T 



A E (x + ie) = e 27re e- 27r »M +)£ (A 1 , A 3 ) 



(4.6) 



M + , e (A 1 ,A 3 ):- L-IJ" 1 ^ 



-Aie 7 ^ + o(i; 




as e — > +00. In particular, for the LE this implies, 
(4.7) L(P,A E )=2ne + L(fi,M + ,) . 



(4.8) 
(4.9) 
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-1 -Aie^ N 
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M + , e (x) = + MW(x), :- l Age _ M ... 

as e — > +oo. Thus, by continuity of the LE, 

L(J3,M +te ) = L(J3,M®) + o(l) 

as e — > +oo. 







In summary, (4.7) and quantization of acceleration yield 



(4.10) L{/3, Af) = 2vre + L(J3, M (0) ) , all e > sufficiently large. 

Since the matrix M in (4.6) is independent of x, the problem is reduced 
to computing the Lyapunov exponent of the constant cocycle (/3,M^). 
Solving the eigenvalue problem for M^°\ we obtain as claimed 

1 + y/l - 4AiA 3 



(4.11) 

(4.12) 
(4.13) 



L(/3,M (0) ) =log 



Using above line of argumentation, the statement about L((3,Af) for e 
large and negative readily follows since then 



A E (x + ie) 



-2ire n 2irix 



M_ 



-l + o(l) -Age - ^ + o(i; 
Aie wi/3 + o(l) 

and L(/3,M(°)) = L(/3,Mi 0) ), where 



(4.14) 



M 



(o) 



(4.15) 



-1 -A 3 e-^ N 
Aie^ 

Case 2 - one of Ai, A3 = 0: Without loss, we assume A3 = 0. 
Taking e — > ±00 we obtain uniformly for x G T, 

-l + o(l) -Aie m/3 + o(l) N 
o(l) 

-l + o(l) (iy 

Aie 7 ^ + o(l) 

The case when e — > —00 being similar, we focus on e — > +00. To this 
end write 



A h (x + ie) =e 2me e 



2irie r , — 2nix 



A E (x + ie) = e -^ e 2vix 



, e — >■ +00 



e —7- —00 . 



(4.16) 



^(x + ^e) =: e^e-^iUx) . 



(4.17) 



Applying Theorem 1.2, we conclude L(f3,R € ) = for e sufficiently large 
and positive. We note that e — > +00 in R e , yields a constant cocycle with 
identically vanishing determinant. 

Making use of (4.16), by quantization of acceleration one concludes 

L(/3, Af) = 2ire , if e > sufficiently large. 
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□ 



We amend that for any choice of A (i.e. in both of the above cases) the Lyapunov 
exponent of (/3, Af) can asymptotically be written in terms of the matrix M, 

(4.18) L(0, Af) = L(/3, M) + 2ix\e\ . 

For later purposes, we note that the matrix M does not depend on A2. 

5. Extrapolation 



In Proposition 4.1 we proved an asymptotic formula for L((3, Af) valid for l arge 
|e|. Using convexity of L(/3,Af) in e, quantization of acceleration and (2.8), we 
obtain the following lower bound for L((3, B E ), 

(5.1) L(/3, B E ) > max{L(/3, M) - /(A), 0} . 

Here we also made use of L(f3,B E ) > 0. We emphasize that in (5.1) we do not 
require E to be in the spectrum. 



For E E E, Theorem 2.10 will allow us to impro ve on (5.1 ) and actually extrap- 
olate L(/3, Bf) to e = 0. Returning to Proposition 4.1, the asymptotic expression 
for L(/3,Af) implies that for |e| > e the acceleration u((3,A E ;e) is ±1. 



Let < e[ < e such that the Theorems 2.5 and 2.10 apply on [— e' 1( e'J. Quan 



tization of acceleration and convexity of L(/3, Af) require u>(/3, Af) to increase on 
-,[] with possible values in { — 1, 0, 1}. 



Moreover, for E E E with L((3,B E ) > 0, by Theorem 2.10 the acceleration 
is forced to jump when passing through e = 0. This amounts to three possible 
situations, a jump in u((3,Af) from — > 1, — 1 — >■ 0, or —1 — > 1. Continuity of 
L(/3, Af) in e, rules out the former two cases resulting in L(f3,Af) = L(f3,M) + 
2n\e\. In particular, (5.1) holds as equality for s uch E. 

/(A). 



0, (5.1) yields L(/3, A 1 



Again 



If, on the other hand, L((3, Bf) 
convexity in e of L(/3,Af) eliminates all possibilities except for 

(5.2) L(/3,Af) = max{/(A),L(/3,M) + 2vr|e|} . 

Since L(/3, M) > /(A) whenever E E E with L(f3, Bf) > 0, we conclude that the 
formula in (5.2) is in fact valid for all A. 



We thus obtain the following expression for the Lyapunov exponent of the 
complexified extended Harper's model, 

Theorem 5.1. For irrational ft, 

(5.3) L(fi,B?) = L(p,A?)-I e {\) , 

(5.4) L(/3,Af)>max{I(X),L(l3,M) + 2n\e\} , 



where equality holds for E E E. L(f3,M) and I e (\) are given in (4-11) and (3.3), 
respectively. 



Finally, taking e = in Theorem |5.1[ Theorem |1.1| is obtained using the follow- 
ing fact, verified by a straightforward computation. 
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Observation 5.1. Denote A := L(f3,M) - I(X). Then, 
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(5.5) 



A 



log 
log 

log 



1 + VI- 


4AiA 3 


2Ai 


1 + vi- 


4AiA 3 



2A3 

1 + V1-4A1A3 



A 2 + v 7 ^ 



I - 4AiA 3 



t/Aj > A 3 , A 2 < A 3 + Ai , 
«/A 3 > Ai, A 2 < A 3 + Ai , 

, if A 2 > A 3 + Ai . 



A > is positive in the interior of region I and A < for A £ II U JJ/. 
Moreover, 

(i) For Ai 7^ A3, A = if and only if A lies on the line segments {Ai + A3 
1,A 2 < 1} or {Ai + A 3 <1,A 2 = 1}. 

(ii) For Ai = A 3 , A = if and only if A G I// or {2A X < 1, A 2 = 1}. 



As a corollary to Theorem |5.1| we obtain a nec essar y and sufficient condition 
for regularity of the analytic cocycle (f3,Af). By (5.2), (f3,Af) is regular if and 
only if A < 0. 

For Schrodinger cocycles, non-regular behavior for zero Lyapunov exponent 
proved to be interesting. In [2J, Avila called such cocycles critical and analyzed 
(lack of) their appearance for generic almost periodic one frequency Schrodinger 
operators [21 H]. 

Motivation for these considerations was the almost Mathieu operator. There, 
the critical point corresponds to the square lattice, separating sub-critical (zero 
Lyapunov exponent and regular) from super- critical (positive Lyapunov exponent) 
behavior (see Eq. (4.4)). 

Implications for spectral theory for general analytic Schrodinger cocyles are 
conjectured (almost reducibility conjecture, see Sec. [7]). 

In analogy to the Schrodinger case we make following defintion: 



Definition 5.2. Consider a Jacobi matrix of the form (1.1) with functions c, v 
analytic and c ^ 0. For given (3, we call the associated analytic cocycle (/3, Af) 
critical if the Lyapunov exponent L(/3, B E ) = and (/3, Af) shows non regular 
behavior. 

For extended Harper's equation Theorem 5^ and Observation 5. 1| yield 

Corollary 5.1. For extended Harper's equation with irrational frequency (3, crit- 
ical behavior occurs 

(i) for Ax 7^ A 3 on the line segments {Ai + A 3 = 1, A 2 < 1} and {Ai + A 3 < 
1,A 2 = 1}, 

(ii) for Ai = A3, in region III and along {2Ai < 1, A 2 = 1}. 
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6. An alternative proof of Theorem 11.11 



As pointed out in the end of Sec. |4j asymptotically L(/3, Af) is independent of 
A2. This interesting fact inspires an alternative proof for Theorem |1 . 1 1 given here. 
First note that above mentioned symmetry is already implied by (4.6) and (4.10) 
without any further computations . 



The arguments of Sec. [5] leading to Theorem |5.1| are independent of the remain- 
der of Sec. J1J this yields a formal expression for the Lyapunov exponent given in 
d5T3b and O). 



We will focus here on the most interesting region III. As we will argue, The- 
orem 

on A2 persists in the limit e — > 0. 



1.1 rests on the following key Lemma which shows that above mentioned 
asymptotic independence of L(fl,B E ) 



Lemma 6.1. Let (3 be irrational. Then, on the spectrum the Lyapunov exponent 
in region III is independent of the coupling A2. 



Proof. First note that for Ai + A 3 > 1, (3.3) implies /(A) = max {log|Ai|, log|A 3 |}. 



Since the matrix M is independent of A2, within region III Proposition 5.1 implies 
that for fixed Ai and A3, L((3,B E ) is constant w.r.t. A2. In particular, in entire 
region III, L(j3, B E ) is determined by its value on the plane {A1+A3 = A2 > 1}. □ 



To see that Lemma 6.1| already i mplie s Theorem |1.1 notice the following: As 
an immediate consequence of Lemma 6.1 we obtain constancy of L(/3, B E ) on the 
spectrum along lines where A3 and Ai are fixed (see Figure [T]). Hence, in region 
III, L(j3, B E ) on the spectrum is determined by the limit when approaching the 
plane Ai + A3 = A2 along such lines from within region III. 

Using continuity of the spectrum in Hausdorff metric as well as joint continuity 



of the Lyapunov exponent in (E, A) (see Theorem 3.1), this limit can equally be 
determined by approaching the plane Ai + A3 = A2 from within region II. 

Based on a duality argument, in [21] we showed zero LE on the spectrum in 
region II for Diophantine /3. For the readers convenience, we give a shortened, 
alternative proof of this fact in Appendix |FJ 

In either case, the duality argument is based on [16] which proves that the 



spectrum in the interior of region I is purely point with exponentially localized 
eigenfunctions. This however is established under the condition that /3 is Diophan- 
tine. As far as the LE is concerned, we can remove the Diophantine condition 



using continuity in the frequency as asserted by Theorem 1.1 



Theorem 6.2. Fix (3 irrational, 
zero on the spectrum. 



Then in region III, the Lyapunov exponent is 



Assembling the pieces, we finally obtain Theorem |1.1[ The Theorems F.l and 
6.2| establish zero Lyapunov exponent in region II and III. Let A G I. Invariance 



Using Theorem 3.1 and continuity of the spectrum in the Hausdorff metric it suffices to 
establish Theorem f.l for Ai, A3 both non-zero. 
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of the density of states under duality (expressed in terms of the map a, see (F.10) 
in Appendix [f]) and the Thouless formula for Jacobi operators imply |20j . 



L(B 



log|c A (a;)|dx+ / \og\E - E'\dn a{x) (E') 



(6.1) 



1 log 




c\{x) 





XI 1 E 



dx + L{B% x f,p) . 



Here, c\ denotes the function c defined in (1.2), whereas in c CT (X) we apply the 
map a to the coupling constants in c. Evaluating the integral in (6.1) for A in 
region I (see Sec. 5 in [20J) yields the formula for the Lyapunov exponent given 
in Theorem II. 11 

7. Implications for general quasi-periodic, analytic one frequency 

Jacobi matrices 

The example of extended Harper's equation not only illustrates a practicable 
method to obtain the Lyapunov exponent (or at least lower bounds) for analytic 
cocycles derived from a one-frequency almost periodic Jacobi matrix but also pro- 
vides useful insights in how the theory originally developed in [2] for Schrodinger 
operators generalizes to the Jacobi case. 

In view of future applications to other Jacobi operators, in the present section 
we shall point out differences and peculiarities of the Jacobi case which relate to 
possible zeros of the determinant of the transfer matrix A E associated with the 
Jacobi operator. 

If not mentioned explicitly, in the following we consider a Jacobi operator on 
Z 2 (Z) of the form (1.1) with functions c,v analytic on Tg, c ^ 0. 

With a give n Jac obi matrix one can associate the two transfer matrices A E and 
B E defined in (2.1). The former gives rise to an analytic cocyle whereas the latter 
is defined only a.e. on T$ due to possible zeros of c(z). In view of the relation 



between the associated Lyapunov exponents, Eq. (2.8), we recall Lemma 2.6 

(1) It is crucial for the extrapolation process discussed in Sec. [5] to understand 
the analytic properties of L(f3,A E ) w.r.t. e. In particular, by Proposition 
|A.l this function is convex in e. Convexity however fails in general for the 
Lyapunov exponent associated with the cocycle ((3,B E ). 

This can be seen on the example of Extended Harper's equation when 
Ai = 0. Then for A 3 = 1 and < A 2 < 1, 



(7.1) 




By Lemma 2.6 and the analytic properties of L({3,Af), L((3, Bf) is still 
piecewise affine, however, since L((3, Af) = 2ir\e\ (Theorem 5.1), convexity 
fails at e = ^MogA2: the right derivative changes from -1 to -2 (with 
increasing e). It can be checked that at this value of e a zero of c occurs. 
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We mention that in general by Lemma 2.6 convexity of L(f3, Bf) may only 
fail at those e where c exhibits a zero. 



In Theorem 2.9 we assumed the analytic cocycle (/3,D) such that det(D) 
is bounded away from zero. That this assumption is indeed necessary may 
be seen when considering ((3, Af) for Extended Harper's equation. To this 
end recall that L(/3, Bf ) = L(f3, (Af)'). For L(/3, Bf ) > 0, regularity of 
((3, (Af)') does in general not imply uniform hyperbolicity if c(z) has real 
zeros: 

Consider Ai = A 3 and 1 > 2Ai > A 2 . Then J e (A) = lo g A 3 + 2ir\e\ in a 
neighborhood of e = 0. Thus, by (4.6) and Theorem 5.1, uj((3, B E ; e) = 



locally about e = even though L((3, Bf) > in the interior of region I. 
Analysis of the Lyapunov exponent serves as an important ingredient in 
spectral analysis. For almost periodic Schrodinger operators with analytic 
potentials the almost reducibility conjecture (ARC) [21 El El U] claims that 
sub-critical behavior of the Schrodinger cocycle on the spectrum implies 
purely absolutely continuous spectrum. 

Section [2] and above remarks suggest to formulate the ARC for Jacobi 
cocycles in terms of the analytic cocycle ((3,Af) (and not for (f3,Bf)). 
However, even for the analytic cocycle (/3, Af) the ARC is false in general. 
For Extended Harper's equation with Ai = and \ 3 = \ 2 > 1, L((3, Af) 



is regular (Corollary 5.1) even though the spectrum has no ac compo- 



nent. The latter follows from Observation 3^ and the following simple 
Proposition. 

We have learnt only recently that this result originally goes back to [11] 
Since the argument is short we include a sketch below. For details we 
refer to [llj. 

Proposition 7.1. Let A such that c(z) has a real root. Then the ac spec- 
trum (T, ac ) is empty. 

Proof. We reduce to the following toy problem: Let 

(Hlp) n := V n 1p n + C n 1p n+1 + C^TlVn-l , 

a Jacobi operator on / 2 (Z) and suppose (c„) ne N is zero infinitely often. 
Then, the Jacobi matrix decouples into finite dimensional blocks whence 
the spectrum is only pp. 

Suppose now c(z) has a real zero, say at xq G T and fix an irrational (3. 
By ergodicity, for a.e. phase 9 the rotational trajectory comes arbitrarily 
close to Xq. Thus for such 6 the operator Hx t e-,p can be considered a trace 
class perturbation of above block Jacobi matrix. Invariance of the ac 
spectrum under trace class perturbations yields the claim. □ 



are grateful to Barry Simon for this reference. 
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Using Observation |3.1[ we obtain a glimpse on the spectral theory for 



Extended Harper's equation: 

Corollary 7.1. Let (5 be irrational. Then, 

(i) For A 3 7^ A 1; S ac = whenever Ai + A 3 = A 2 . 

(ii) For A3 = Ai, S ac = along 2A3 > A 2 . 

Above example shows that, in general, ARC fails if c(z) has zeros. We 
do however believe that the conjecture is still true whenever detAf(x) is 
bounded away from zero: 

Conjecture 7.1 (ARC for Jacobi cocycles). For an analytic almost peri- 
odic Jacobi cocycles with determinant bounded away from zero, sub-critical 
behavior of (f3,Af) implies (purely) absolutely continuous spectrum. 



Appendix A. 
A.l. Convexity of L(0, D £ ) w.r.t. e. 

Proposition A.l. Let ((3,D) be an analytic cocycle in the sense of Definition 



2.1, Then, L(/3,D e ) is convex in e. 



Proof. The main point here is to convince the reader of convexity for singular ana- 
lytic cocycles. For D G C W (T, SL(2, C) (or more generally, for D G £ W (T, M 2 (C))), 
it is immediate that L(/3, D € ) is subharmonic w.r.t. e (by e.g. a Craig-Simon type 
argument [8]), in which case convexity easily follows, as pointed out in [2j. We 
mention that the argument presented here can also be generalized to analytic 
cocycles on T d , d > 1. 

Since convexity is preserved under pointwise limits, by the definition of the LE, 
(2.5), it suffices to verify that 

(A.l) - I log||D (n) (a; + ie)||djc , 

n Jj 

is a convex function in e for each n G N. 

Considering M 2 (C) as equipped with Hilbert-Schmidt norm, this readily follows 
from: 

Lemma A.l. Let {fj}f =1 be a finite sequence of functions holomorphic in some 

neighborhood of S 1 , and set := log J2f=i\fj\ 2 - Then, J(e) := j' T ,/i 0O (e 27ria: e -27re )ds 
is convex in e. 

Proof. Follows by a simple approximation argument: For 77 > 0, let 

h v := log ^X^il/jl 2 + Vj an d correspondingly, J v (e) := f T h v (e 2mx e~ 2mt )dx. 

Clearly, J v \ Joo, as ^ -> + , whence the claim is reduced to J v . 
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Direct computation shows that h n is a smooth, subharmonic function in a neigh- 
borhood of S 1 . E 

Finally, considering e as a complex variable, the identity 

(A.4) J?l = 47r 2 e-^ Ree ^ 



dede dzdz 
implies that J v is a smooth, subharmonic function in e which is constant in Ime, 
i.e. convex in Re e. □ 

Remark A. 2. Since subharmonicity is preserved under monotone decreasing limits, 
the proof of Lemma A.l also shows that J v {e) and hence L(/3, D e ) is subharmonic 
in e. 

□ 

Appendix B. Proof of Lemma 12.31 

Stepl: We first establish the statement for trigonometric polynomials, where 
g can be computed explicitly. Since, 

n 2n 

(B.l) cje 2wijx = e~ 2ninx Cj- n e 2mjx , 

j=-n j=0 

it suffices to consider trigonometric polynomials which are monic algebraic 
polynomials in the variable e 2mx of degree r > 1. In particular, / can be 
factorized according to its roots (counting multiplicity), 

r 

(B.2) / = Y[ (e 2nix - e 2 ™^+^)) , 

j=i 

where for some 6 > 0, Xj G [0, 1) and \ej \ > 5 > 0, 1 < j < r. 

Since / is holomorphic and bounded from zero on the strip |Imz| < 5, 
there exists a holomorphic function q defined on |Im,z| < 5 such that 

(B.3) e* = f ; 

in particular, q satisfies 

(B.4) q'(z) = f j^ ,\lm(z)\<5 . 



If {<7j} is a finite sequence of holomorphic functions with ~^2j\gj(z)\ 2 ^ 0, one computes: 

j 2^j=i\9j\ y^k li<k 

(A.3) > 1 J2(\9iW9k\-\9'k\\9i\f>0- 
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On the other hand, for i£l, 



<b*) m-t 



(B.6) 



Considering the individual summands in (B.5), we compute for 1 < j ' < 
r and 16K 

2me 27zix ( 27r iJ2™ =0 e 27Tix J n e- 27Tne >e- 27Tinx if ej > , 



g2-7r«x Q2iri(xj+iej) j 27T2 e~ 27rtTKB J -g27rejng27rinx jf e . < Q 

Hence, for x G IR, q is given by, 



fit) 

(B.7) - g(0) = / ^yfyd^ = 27rm + x + h(x) 

Jo J [t) 



where n + represents the number of j G {1, . . . ,n} with ej > 0, and h G 
Cg (R/Z; C). Notice by construction, g(0) = log/(0), where the branch of 



the log is fixed by (B.3). 

Let now g(z) := ea 9 ' 2 ', which by construction is holomorphic on |Im(z)| < 



5. For iGl, the explicit expression in (B.7) already implies that g(x) is 
2-periodic satisfying g 2 (x) = f(x). Upon use of the uniqueness theorem for 
holomorphic functions, both these properties extend to all of |Im(,2)| < 5. 
Step 2: If / G Cf{R/Z; C) is not a trigonometric polynomial, let 5 > such 
that f(z) 7^ 0, \/z G TV Approximate / uniformly on T$ by trigonometric 
polynomials p n . 

From step 1, we obtain holomorphic q n on |Im(z)| < 5 such that 

(B.8) g n := e^ G C?(R/2Z; C) , g 2 n = p n and 

(B.9) g n (0) = logp„(0)^log/(0) , 

defined (eventually) with respect to a common branch of the log since 
/(0) ^ and p n -> /. 

Since g n satisfies the differential equation = — on |Im(2;)| < 5, and 
by construction p n — >• / uniformly with / bounded from zero on TV, (q' n ) 



is easily seen to be uniformly Cauchy. Using (B.9), this in turn implies 
uniform convergence of (q n ) on T§. 

Finally, letting q := lim n _>.oo9n an d defining g := e^ q we obtain the 
lemma's claim. 

Appendix C. Proof of Lemma 12.121 

The proof for the cocycle (A E ) ; ) will at the same time prove the statement 
for ((3,B E ). Again, the function log|c(z)| is harmonic on the strip |Im,z| < 5, 
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some 5 > 0. Harmonicity and Birkhoff's ergodic theorem imply that Vx 6 T 

->0 , 



(C.l) -log 



?? 



3=0 



|det + i/3)|V2 



as n — >■ oo. 

Suppose (/3, (A )') is uniformly hyperbolic for E G S. According to Definition 
2.8, let be the unstable/stable directions for the positive half-line 

solutions of H x ip = Eip, iGl Note that uniform hyperbolicity on the positive 
half line implies uniform hyperbolicity on the negative half-line with u + = and 

M_ = S + . 

Let w := (ip(0), tp(— 1)) G C 2 , ||w|| = 1, be an arbitrary initial condition. For 
x G T arbitrary, choose unit vectors u x G u + (x) and s x G u_(x). Then, if = 
+ a x s x for some complex a x , (3 X . Without loss of generality, we may assume 
Px 7^ (otherwise apply the argument for iteration in the negative direction, since 
then v = a x s x G U-(x)). 



Using flail ), 

o 

log|| ]"[ (A B )'(x+i/3)|| = -log|| n B E (x+jP)\\+o(l) ,asn^oo. 



1 ° 1 



j=n—l j=n—l 

Therefore, uniform hyperbolicity of ((3, (A )') implies that the solution of H x ijj = 
Eip associated with the initial condition w increases exponentially in at least one 
direction, w and x were chosen arbitrarily, hence Vx G T, E is not a generalized 
eigenvalue of H x . Since uniform hyperbolicity is an open condition, using above 
argument again, Vx G T, E cannot be a limit point of generalized eigenvalues of 
H x either. Hence, we obtain E is not in the spectrum, a contradiction. 

Appendix D. Upper-semicontinuity of the LE; Proof of Theorem 
Since all matrix norms are equivalent, following it is convenient to equip M2(C) 



with the Hilbert-Schmidt norm. Note that Remark 2.16 also applies to log + and 
continuous M 2 (C)-cocycles, i.e. for every n G N , (f3,D) h-> J T log + \\D^ n \x) \\dx 
is continuous in C° (T, M 2 (C) ) . 

On the other hand, by Fatou's lemma, (/3,D) >-)■ Llog_||Z)W(a;)||da; is upper- 



semicontinuous. Recalling (2.5), the statement of the theorem follows since the 



infimum over upper-semicontinuous functions is upper-semicontinuous. 

Appendix E. Proof of Lemma 12.171 

For K C (0, 1) compact, let m := minf^ and M := max.fr. Fix / G J. By 
basic properties of convex functions, 

(E.1) D +(9)( M) < ^< + »)-^) , 
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for any g E J 7 and < h < 1 — M. 
In particular, 

(E.2) inf SUP D Mi M)< fiM + h) - f<M) , 

* >O ll0-/ll[o,i]<* h 

for any < h < 1 — M. Thus, taking h — > 0+ we obtain: 
(E.3) inf sup D + (g){M) < D + (f){M) . 

5>0 ll9-/ll[o,i]«5 

A similar argument shows, 
(E.4) £)_(/) (m) < sup inf D-(g)(m) . 

<5>0 ll3-/ll[o,i] <<5 

Recalling that D±(g) is monotone increasing on [m, M], we obtain the claim. 

Appendix F. Duality 



In view of the alternative proof of Theorem 1.1, we establish the missing link 



in Theorem 16.21 

Theorem F.l. Throughout region II, L((3, B E ) = for all E in the spectrum and 
irrational (5. 



Since the Hamiltonian in (1.1) generalizes the almost Mathieu operator, it 
should not come as a surprise that its spectrum expressed in terms of the density 
of states exhibits similar symmetry w.r.t. certain permutations of the coupling 
constants as the almost Mathieu operator. In fact, this so called duality has been 
explored for the present model in previous works [201 HS1 Ell EEH ES]- A more 
general discussion of duality can be found in [20] . 



In this Section we will prove Theorem F.l, adapting an idea of Delyon [TO] , 
that originally was used to establish absence of point spectrum in the sub-critical 
region of the almost Mathieu equation. 

Let a denote the map, a(X) := ^(^3, 1, Ai), where, as earlier, A = (Ai, A2, A3). 

Referring to the partitioning of the parameter space introduced in Sec. [TJ the 
map a acts as follows 

(F.l) o{I) = II , a(III) = III . 

In this sense, regions I and II are dual, whereas region III is self-dual. 

In [16] (see Theorem 1 therein) it was proven that for Diophantine (5 and A in 
region I, the spectrum of is only pure point (i.e. for a certain set To of phases 
9 with I To I = 1). Fix A in region I and ^ G To. Let <pe be a normalized eigenvector 
of H\ e with associated eigenvalue Eg. If we let 



(F.2) Mx) = J2Mn)e 



2irinx 
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-.2m9m 



<Pe(v + /3m) \ defines a random sequence indexed by points in T. Set 

(F.3) MV, m) := e 2mdm M V + (3m) . 

Since for each fixed m £ Z, ipe(.,m) £ L 2 (T), direct computation shows that 



c<t(a)(- + + m + 1) + c (T ( A) (. + /3(m - 1))^(., m - 1) 

(F.4) +u(. + £m)V>e(.,m) = X^Ee^^m) , 

for each m £ Z. 

Equation (F.4) holds in the sense of L 2 -functions, whence it is true on a set 
Q m of full measure in T. In particular, letting Q := n me zfi m , for each T] £ Q the 
sequence {i/)o(rj, m)} meZ satisfies the finite difference equation 

(F.5) H a{x)iV ip e (m .) = X^Eeipeiv, •) • 

At the same time, for e > normalization of <pg yields 

2 



(F.6) 



n 2 



< OO 



Thus for every e > and Lebesgue a.e. r\ there exists C e (rj) > such that 
(F.7) \Mv,n)\ < C e (?7)n 1/2(1+e) , Vn £ Z . 

It is a general fact for a family of bounded ergodic Jacobi operators 
(O, J 7 , P) a probability space, that fixing an energy i£ [5], 

(F.8) P{£ is eigenvalue for # w } = . 

Ma king use of a version of Osceledets Theorem for Jacobi operators, Theorem 
G, (F.8) implies that for P-a.e. u any ( £ C z satisfying H^C = E( has to 



exponentially increase in one direc tion with rate given by the Lyapunov exponent. 

' , (3) =0. Hence, if 



The sub-exponential bound in (F.7) then implies L IBJ 2 ^ 
we let 



U M# 

6>eT 



A;0j 



^(A)^ 



for all E £ 21. In fact, continuity of the Lyapunov exponent 
w.r.t. the energy yields 



(F.9) 



LiB;if,/3 



, ME £ 21 . 



Set 93 := R \ 21. 

Denote by dn ({^,6*} ; E) the density of states for the family of ergodic oper- 
ators {H\ } o}. As shown in a more general context in [20], duality preserves the 
density of states. For the present model this yields, 

(F.10) dn {{H Xt0 } ; E) = dn {{\ 2 H a{xl g} ; E) . 
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Note that the spectrum of {H\ t g} is contained in 21, hence n ({H\ e } ; 03) = 0. 
Using (F.10), we obtain 

(F.ll) Z({\ 2 H a(x) }) C2T. 

Finally, the spectral mapping theorem, (H, and (|F.ll|) imply L(5f (A) ,/3) = 

for any E G S ({A 2 i^ (7 (A)})- Since A was chosen arbitrarily in region I, this implies 
Theorem F.l| for Diophantine /3. 

Finally, as mentioned in Sec. |6l the Diophantine condition can be removed 



using continuity in the frequency as asserted by Theorem 3.1 



Appendix G. Oseledets Theorem for Jacobi operators 

In Appendix [F] we made use of Oseledets theorem to relate the Lyapunov expo- 
nent to the exponential decay rate of solution to the Schrodinger equation. The 
authors noticed that in the literature for Schrodinger operators, Oseledets the- 
orem is usually given as a theorem for SL(2,C) matrices even though Oseledets 
original statement allows GL(2,C) [22] . 

The relation of one to the other is, of course, just a matter of appropriate re- 
normalization of the determinant. We include the simple argument below and 
show its consequences for ergodic Jacobi operators. A more general, deterministic 
statement which even allows for not invertible matrices can be found in 



Theorem G.l (deterministic formulation). Let (A n )neN be a sequence in GL(2,C). 
Setd n : = (detA n )a, D n := j^A n , n G N. Suppose thatm := lim^oo ^ XlLi lo S K 
exists and is finite, lim^oo - log ||Ai|| =0, and that L = L(A) := lim^oo i log \\A n 
exists and is finite with L — m > 0. Then, there is a one- dimensional subspace 
SCC 2 such that for v G C 2 



(i) ilog||A n ...A 1 u| 

(ii) i log ||A n ... AH 



-> —L + 2m, ifveS, 
\L, ifveC 2 \S. 



0. 



Proof. Since m exists and is finite, - log \d n \ - °°> 0, so lim^oo ^ log ||-D n || 
Note that L — m = lim^oc - log \\D n . . . Di\\, hence by the SL(2,C) version of 
Oseledec theorem, there exists a one-dimensional subspace SCC 2 such that for 
v G C 2 



n 



1 1 n 

log \\A n . . . Aiv\\ = - log \\D n . . . Div\\ + lim - ^log \d, 
n n^oo n — ' 



k=l 



-> — L + 2m , for v G S 



1 1 1 n 

-log||A n ...^|| = -log\\D n ...D 1 v\\ + lim - Vlog|tZ, 
n n n^-oo n — ' 



k=l 



% L , for v G C 2 \ S 



□ 
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Remark G.2. (i) Let A : Q — > GL(2, C), where (Q, J 7 , P) is a probability space 
and T : Q Q ergodic. If A n (x) = A(T n x), n e N, then by Birkhoff's 
ergodic theorem the Ceasaro mean m becomes, m = E(log|c?i(a;)|). For an 
ergodic setup, we may hence replace the condition on finiteness of m by 
E(log|di(x)|) < oo. 
(ii) Consider a bounded ergodic Jacobi operator, i.e. 

(H^) n := b{T n u)i) n + a{T n u)^ n+1 + a(T n_1 w)-^n-i , 

where T is ergodic on the probability space (fl, J 7 , P), a : Q — >■ C \ {0}, b : 
Q — > R are bounded random variables and E w log|a(w)| < oo. 

In this case, above sequence {A n ) n£ i is the associated sequence of one- 
step transfer matrices, 

1__ (b{T n u) - E -a{T n ~ l uj) 

a(T n u) V <T n u) 

where BeR. Then, since T is measure preserving, (i) implies m = 0. 
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